General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



N76-28738 


(NASA-CR-1 48533) SOLAR MICROCLIMATOLOGY 
Final Report (Helio Associates, Inc., 

Tucson, Ariz.) 71 p HC $4.50 CSCL 04B 

Oncl as 

G3/47 15184 


SOLAR MI CROCLI NATO LOGY 

CONTRACT NO, NASW-2745 


FINAL REPORT 
NOVEMBER 1975 


Prepared For 

NASA Headquarters 
Washington, D. C. 20546 


Prepared By: 


Dean B. McKenney 
William T. Beauchamp 




H2LI0 ASSOCIATES, INC. 

P.O. Box 17960 
Tucson, Arizona 85731 


' BIBLIOGRAPHIC DATA i* Report No. 
SHEET 


dw Tftlc and Subiuic 


"Solar Microclimatology" 


ATTAr^ lT..._J 

3. Recipient’s Accession No. 

5. Report Date 

November 1975 


7. Author(s) 

W. T. Beauchamp and D. B. McKenney 


9, Perform ing Organization Name and Address 

/ 

Helio Associates, Inc. 
P. 0. Box 17960 
* Tucson, Arizona 85731 


12. Sponsoring Organization Name and Address ’ 

RASA Headquarters 
Washington, C. D. 20546 


8. Performing Organization rtept. 
No. 


10. Projcct/Task/U'ork Unit No. 


11. Contract/Grant No. 

NASA 2745 


13. Type of Report £: Period 
Covered 

Final Report Nov 1975 


14, 



16. Abstracts 


Insolation data was collected and analyzed for a network of 
sensors distributed over an area of several square kilometers. 

The analyses of this data yielded probability distributions of 
cloud size, velocity, and direction of motion which were compared 
with data obtained from the National Weather Service. Microcli- 
matological analyses were also performed for suitable mpdeling 
parameters pertinent to large scale power plant design. • ■' 


17, Key Words and Document Analysis. 17a. Descriptors 

Solar Energy 
Solar Insolation 
Solar Power Plant Design 
Cloud Motion 

Insolation Instrumentation 
Insolation Measurements 


17b. Identifiers /Open-Ended Terms 


Page rq 
op P °0R Q 


17c.- COSATI Field/Group 


IB. Avnilabilicy Statement 


Release Unlimited v 


19*. Security Class (Ti:is 
Report) 

20. Security Class llliii 
Pace 

itKic’i Kcimcti 


21. No. of "aces 


I 22 , Price 








ABSTRACT 


Insolation data was collected and analyzed for a network 
of sensors distributed over an area of several square kilometers. 
The analyses af yielded probability distributions of 
cloud >sSsg£2££on of motion which were compared 

with data ob tain^d^rrg ^ jgfhe National Weather Service. Micro- 
climatological analyses were also performed for suitable modeling 
parameters pertinent to large scale power plant design. 


fX? 141 ' PAGE is 
f poor quau^ 




TABLE OF CONTENTS 


Abstract , ± 

SECTION 1.0 INTRODUCTION i 

SECTION 2,0 FIELD AND INSTRUMENT DESCRIPTION . ... 3 

2.1 Instruments - Physical ......... 3 

2.2 Instruments - Electronics . g 

2.3 Calibration and Installation ...... 3 

SECTION 3.0 DATA ANAYLSIS ............. 12 

'3.1 Event Anaylsis 12 

3.2 Limitations on the Data . .... . . . 21 

3.3 Reduced Data - Velocity, 23' 

Angle, and Size 

3.31 Correlated data . , . . *. , . . . 23 

3.32 Correlated plus 

uncorrelated data . . . . ' . ... 27 

3.4 Power Plant Site Analysis . . .... .28 

3.5 Fourier Representation „ . . . . . . . . 37 

3.6 Calendar Representation . . . , ... . 43 

3.61 Comparisons with NSW data . . . . 45 


SECTION 4.0 CONCLUSIONS . , .49 

SECTION 5. 0 RECOMMENDATIONS FOR CONTINUED RESEARCH . 50 
SECTION 6.0 APPENDICES . . . . . . . ....... .52 

Appendix I - Figures 20 through 27 . . . 53 
Appendix II - References . . % .. .. . . . 62 
Appendix III - Acknowledgements . . . . 53 



FIGURES 


Figure 1. Sketch of the Moving Bar Occulting 

Radiometer (MB OR) 4 

Figure 2. Detector design for the MBOR 5 

Figure 3. Schematic diagram of Helio circuit 

used for the current- to -voltage 
transducer 7 

Figure 4. Output of MBOR Instrument 9 

Figure 5. Location of instruments to measure 

spatial distribution of solar inso- 
lation 11 

Figure 6 . Analog output of MBOR for a clear day . , .13 

Figure 7 . Analog output of MBOR for a typical 

partly cloudy day 14 

Figure 8. Schematic for the measurement of the 

time shift for a single event at two 
stations 16 


Figure 9. Schematic and modeling for the velocity 

calculations for the field of 5 detectors . 17 


Figure 10. An example of the components of the 

measured velocity (V r ) due to the actual 
motion of the cloud (V m ) and due to the 
shape of the leading edge of the cloud 
(V f ) 20 

Figure 11. Velocity distribution determined from 

spatially correlated data 24 

Figure 12. Direction rose for correlated data 25 

Figure 13. Distribution of calculated cloud sizes, 

D, (kilometers) from spatially correlated 
data 26 

Figure 14. Velocity, (speed) distribution for cloud 

motion over the field using all of the 
data 29 


iii 



30 


Figure 15. Direction rose based on all of 
the data 


Figure 16. Frequency vs the time any cloud 
is over an elemental area, AA, 
in the field 33 

Figure 17. Frequency vs the percentage of the 
field area that is occulted by an 
individual cloud v .... 35 

Figure 18. Frequency vs down time for passage 
of cloud over whole field - largest 
dimension 36 

Figure 19. The Fourier spectra for clear and 

partly cloudy days 39 

Figure 20. Calendar display of the analog data 

for the month of November 1974 54 

Figure 21. Calendar display of the analog data 

for the month of December 1974 55 

Figure 22. Calendar display of the analog data 

for the month of January 1975 56 

Figure 23. Calendar display of the analog data 

for the month of February 1975 . . . . . .57 

Figure 24. Calendar display of the analog data 

for the month of March 1975 . 58 

Figure 25. Calendar display of the analog data 

for the month of April 1975 59 

Figure 26. Calendar display of the analog data 

for the month of May 1975 60 

Figure 27, Calendar display of the analog data 

for the Month of June 1975 ........ 61 

Figure 28. Percentages of days elassifed as clear, 
partly cloudy, or overcast for the 
months of November 1974 through dune 
1975 (Cloud Cover) 46 


Figure 29. Percentages of days classified as clear, 
partly cloudy, or overcast for the months 
of November 1974 through June 1975 
(Possible Sunshine). ........... 47 


iv 



Table 1 

Table 2 
Table 3 

Table 4 


TABLES 


Table for assignment of compass 
' calculations from the angle of 
motion determined in Fig. 5. . , . •* ‘ , 

Microclimate lo gy questions pertinent 
to power plant site analyses . . . r . 

Day classifications and the results 

of Fourier plus ■ phase modeling for 

each ' . .■ ‘ . ", . , . . . . . 1.1 . I: ■ . . . ' 

Division of Calendar-Day data into 
three classifications showing simi- 
larity of the eight Winter/ Spring/ 

Summer months * 



SECTION 1 . 0 INTRODUCTION 
* 

It has become apparent in recent years that Solar Energy 
can be used for electric power production by several methods. 
Because of the diffuse nature of the solar insolation, the area 
involved in any central power plant design would encompass several 
square miles. Detailed design of these large area collection 
systems will require precise knowledge of the local solar insola- 
tion. Information regarding the total available insolation and 
suitable integrals 'of the data are currently available with off 
the shelf instruments and from many National Weather Service in- 
stallations, However, for many system designs, detailed informa- 
tion will also be needed concerning the temporal nature of the 
insolation and the local spatial distribution. Such information 
is not currently available in the literature. Visual observa- 
tions show that it will be highly probable that different portions 
of a large - several km across - collecting area will experience 
different insolation characteristics at different times of the 
day due to clouds smaller than km dimensions. It will be necess- 
ary ; fo take these effects into account in future analyses of solar 
collector systems , 

The purpose of -this research project was the collection 
and analysis of insolation data from a network of sensors distri- 
buted; over an area several kilometers on a side. Both direct- 
normal and total-horizontal insolation data were collected in 
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order to model various solar absorber system types. The field 
of detectors used to gather the data formed a 3.6 x 1.6 km grid. 
This grid allowed for two independent evaluations of the velocity 
and direction of motion for clouds passing overhead. The calcu- 
lated velocity was used with the duration of the occiiltation to 
obtain a typical dimension for each cloud. A Fourier modeling 
of temporal and spatial distribution of the solar insolation over 
the field was also done. A complete discussion of this modeling 
and the microclimate trends are presented in the following sec- 
tions'. A brief summary is presented in the next paragraph. 

The data analysis reported here has yielded much valuable 
information despite the amount of varience in cloud types observed, 
the variety of weather conditions in the area, and the complexity 
of the resulting "events”. Several interesting features are found 
for Tucson’s microclimate. These features readily illustrate 
the information available for large scale system designs. The 
velocities measured for the occulting clouds are principally in 
the range from 0.60 km to about 9 . 7 km with the most probable 
extent in the direction if motion of about 0.19 km. The data span 
six months of operation and include more than 200 individual 
events. About 40% of these events were represented by sufficient 
data to calculate two independent values of the velocity and 
direction of motion. About 67% of these eorrelatable events 
showed good agreement in both values of the velocity and direc- 
tion. This suggests a 67% confidence in any single vector velocity 
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calculation. The analysis results are discussed in detail in 
Section 3. A brief review of the instruments and the field 
layout follows. 

SECTION 2.0 FIELD AND INSTRUMENT DESCRIPTION 

2.1 Instruments - Physical 

The instr um ents used to specify the solar insolation are 
modified pyranometers with a shadow bar that interrupts the direct 
solar beam as it rotates around the detector aperture, A sketch 
of the Moving Bar Occulting Radiometer (MBOR) is shown In Figure 
1, The compact size of this instrument - about 20 x 6.5 x 7 cm - 
makes it suitable for field operation. A 4.8 mm diameter circu- 
lar detector aperture is located on the axis of revolution of a 
circular occulting bar. The surface of the aperture is opal 
glass. A detector is mounted beneath the aperture plate as 
shown' in Figure 2. The area Immediately under the aperture 
plate houses, a neutral density filter required to keep the 
detectors from saturating at the maximum f lux levels . This 
surface- filter- detector combination insures . an output that is a 
linear function of the flux level and a cosine like angular res- 
ponse. The occulting bar is driven at a rate of -J- RPM, and has 
a cross section such that a sharp shadow just- covers the aperture 
once per revolution of the bar in the field of view. To minimize 
the error in the diffuse flux caused by the presence of the bar, 
Its width is specified by 

w = d + 2R sin (%°) 

. -3- .. 




Figure 1. Sketch of the Moving Bar Occulting Radiometer (MB OR) 
a ~ .aperture Plate 

b “ Shade bar •- 

c - Motor housing . 

- d - Detector 






Figure 2. Detector design for the MB OR 

a — Aperture plate 
b - Diffuser 

c ~ Detector locating- sleeve 
d - Cell cylinder 
e •• Detector (TO- 5 can) 
f - Transistor socket 



where d is the aperture diameter and R is the radius of curvature 
of the bar. The bar is also cut so its arc length is sufficient 
for the maximum solar altitude. A correction factor has been 
calculated which takes into account the diffuse radiation lost 
to the detector due to the presence of the bar in the 2^ stera- 
dian solid angle of view. When the bar is out of the detector 
field of view, the unit measures the total flux. While the 
shadow is over the -detector , the instrument is recording the 
diffuse insolation. 

2.2 Instruments - Electronics 

The electronic circuit used to amplify the detector out- 
put for Chart Recorder display is shown in Figure 3. This opera- 
tional amplifier circuit acts like a current to voltage converter 
for the silicon detector operating in a zero bias photovoltaic 
mode. The load resistance for the photodiode is near zero due 
to the characteristics of the circuit. This guarantees a linear 
response for the device in terms of output voltage vs flux level. 
The resistor R. is used for null adjustment, of the amplifier* and 
the gain is adjustable through R x and R 2 . A reference current 
input is also obtainable from V re f through the precision resistor 
R to allow for periodic gain stability calibration checks. This 

Jr 

precision reference current source was external to the five indi- 
vidual output circuits employed in this research, and was used 
as a common source for calibration of the electronics in the field, 

- 6 - 




C ~ >0.47yf R 1 - 2MQ R 2 *= 2MJ2 P 3 = 12M2 


Figure 3. Schematic diagram of Helio circuit used for 
the current-to-voltage transducer. 

Pj , R 2 feedback resistors (R^. = Rj + J?2 or ^1 J * 

J ?3 zero offset adjustment. 

V _ reference voltage for circuit gain calibration, 
ref 

R reference precision resistor . 



Our experience with these detector-circuit combinations over 
the past two years has shown them to be extremely reliable and 
trouble free. In addition, the gain was found to be stable 
within 1 % of full scale and the zero drift was less than 1 % 
per month. 

The chart recorders used for the read out of the data 
are Soltec 10" potentiometric recorders Model 251. The chart 
speed was l"/hr' to * start and 2"/hr in later months. This allowed 
sufficient reading accuracy for measurements of the flux qualities 
at 6 min intervals. Each recorder was equipped with an event 1 
marker driven by an hourly clock pulse derived from a Bulova 
interval timer Model TE240. An example of the resulting output 
record is shown in Figure 4. This analog output shows, at a glance, 
the total, the diffuse, and the direct, § d , insolation 

as indicated. 

2.3 ■ Calibration and Installation 

The instruments were calibrated with respect to an Epply 
Pyrheliometer by measuring the direct normal solar flux. Each 
MBOR was tilted perpendicular to the sun, and readings of the 
total and the diffuse flux were taken. The difference in these 
numbers was then compared to the reading from the pyrheliometer 
and the sensitivity calculated. The pyrhe Home ter was recently 
calibrated against a group of reference instruments traceable to 
the National Bureau of Standards. All five MBOR's were adjusted 
electronically to have the same sensitivity of 7.50 mv/w /^2 • 
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Measurements were also made of the minimum resolvable 
time for the duration of occultations and the difference in time 
of occurance of these events. It was found that 0.1 inch is the 
smallest resolvable occultation due to the width of the chart 
recorder pen. For a 2 M /hr chart speed, this results in a maxi- 
mum resolvable occultation period of 0.3 minutes. As is discussed 
later, "events" of this short a duration were difficult to measure 
because they could not be clearly identified and correlated on all 
five records . The instruments were located on the grid shown in 
Figure 5., In this grid, the Helio laboratory was Station 5. The 
other locations were small convenience markets where electricity 
and some equipment security were readily available. Relative 
timing between events (occultations) are used to specify the 
spatial or velocity profile of clouds over the field. Measure- 
ments of the event timing, between Stations 1 and 2 and between 
Stations 3 and 4, yield two independent measurements of the X 
component of velocity. Once the velocity Is specified, the size 
of the cloud in the direction of motion can be calculated by 
examination of the duration of the individual event at each of 
the five stations. The results of these calculations are values 
for the representative size, speed, and direction of motion of 
the clouds. 
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Figure 5. Location of Instruments to measure spatial 
distribution of solar insolation. The 
asterick represents the Helio Laboratory. 
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SECTION 3.0 DATA ANALYSIS 


The analog insolation data records obtained from eacn 
station in the field of detectors are of the form shown by the 
clear day record in Figure 6 and the cloudy day record in Fig- 
ure 7. These data are presented in the following consolidations 
or reduced forms: 

1. ) Event analysis employing the five MBOR 

instruments ; 

2. ) Calendar day display of the MBOR data 

for one site; 

3. ) Frequency or temporal analysis - Fourier 

Analysis - and modeling for the MBOR data 

as read digitally for selected days. 

Details of these data presentations are discussed in the follow- 
ing paragraphs. 

3 . 1 Event Analysis 

For the purpose of this analysis, an event is defined as 
any sudden change in the insolation level. On a clear day, the 
only changes in the insolation are due to angular effects and 
atmospheric attenuation due to longer path lengths . These 
changes are gradual and result in a smooth data graph such as 
shown In Figure 6. On cloudy days, however, there are sharp 
changes In the insolation as clouds occult the sun or holes in 
the clouds briefly expose the detector to sunlight. The dura- 
tion of these events are typically between 1 and 60 minutes, 
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and they can occur frequently and rapidly throughout the day. 

The duration of each event and the velocity with which the edge 
of the cloud or hole traverses the large collector field are im- 
portant considerations for the design of solar power plants. 

Each individual detector in the field used in this study 
can yield the mean duration of the event and, therefore, the 
duration of the loss of input to a part of the collector field. 
Comparison' of the events between stations can yield the vector 
velocity of motion of the cloud across the field, To evaluate 
these quantities for selected cloudy days, we measured the event 
parameters shown in Figure 7. The quantities T D and r £ u are the 
times for the down transition and the up transition for the event 
respectively. From these numbers a mean difference or shift in 
the time of occultation between adjacent stations is calculated. 
These differences, shorn schematically in Figure 8, are used to 
specify the velocity and direction of motion of the occulting 
cloud as shown in Figure 9. Both types of events, ’’hole" and 
"cloud 15 , are represented in Figure 7. The two voltages, V y and 
V shown were also recorded for internal use.’ If the difference 
between these levels is not a large fraction of the total, then 
the event is probably uninteresting for anlysis because it would 
not be representative of a significant loss in input energy to the 
collector system.. Therefore, we put a limit on the depth of the 
occultation or the height of the peak such that (V r - V B )/V T >0.2. 
This threshold value will preclude consideration of the small 
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Figure 8. Schematic for the measurement of the time 
shift for a single event at two stations. 


ht D Lt V 



At„ ® Tams shift for the start of the occultation . 

D 

Aty -■=• Time shift for the end of the occultation, . 

V^r V , V are the signal voltages at the start , the end and the 
L B bottom of the event respectively. 
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Figure 9 . 


Schematic and modeling for the 
velocity calculations for the 
field of 5 detectors „ 
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The L . . ' s are the distances between stations along the X £ X direction. 
The At. ,'s ere as shown in Figure 4 , 
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variations in insolation such as characterized by the top of 
the relatively smooth portion of the data in Figure 7 around 
11:00. The angle, 6, -will be converted to compass direction 
according to the schedule in Table 1. 

The Velocity vector obtained is an approximation to 
the actual direction of motion due to the possible odd shape of 
any one cloud. As an example of the reason for this approxima- 
tion, consider the hypothetical cloud shown in Figure 10. Be- 
cause the edge of the cloud is tilted to the direction of motion, 
there is an apparent velocity vector normal to the edge, V^. 

This component combines with the actual velocity of motion 
to form what we see as the velocity, V r . The magnitude of the 
difference between V r and Vjj, depends greatly on the shape of 
the particular cloud observed, and although probably small, is 
not easily evaluated. Therefore, the velocities quoted in this 
report will be approximate, but still very, valuable in power 
plant design. 

The event duration can be combined with the velocity data 
to obtain representative cloud size data in the form 

5 ip rri 

g = lj x Uri D,j 

5i=j V 

where V = (V 2 + V 2 ) z . These cloud ’’dimensions" are representa- 
x y 

tive of the extent of the cloud in the direction of motion only. 
This is the dimension of principal interest for loss of insolation 


- 18 - 



Table 1. Table for assignment of compass calculations 

from the angle of motion determined in Fig. 5. 









input considerations. The results of these data analyses are 
presented below. 


3. 2 Limitations on the Data 

Because of the vagarious nature of the weather and cloud 
formation, the analyses used here are necessarily limited and 
approximate. In spite of the complexity of the microclimate 
modeled, we are able to make significant conclusions based on 
what we feel is a good reliability figure for the data. During 
the period of observation, we were able to interpret 225 "events” 
common to at least three stations in the field and subject to the 
restrictions discussed above. Of these, 134 events did not have 
usable data from all five stations, and could not he used to 
calculate the two independent values of velocity. One of the 
principal reasons for the loss of comparable data was the late 
speed change in the home base chart recorder. This recorder 
was maintained for a longer period of time at the 1" per hour 
chart speed convenient to the calendar data display presented 
in a later section of the report. This loss of Station 5 data 
eliminated 32 events from the comparative data.. Ninety -three 
events yielded measurements for all five stations so that the 
two independent velocities could be calculated. Of these 93 
calculations, 62 events yielded agreement in the calculations 
of both the velocity magnitude and the direction of motion. 

For the purposes of this analysis, the difference assumed as 
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limits on the "agreement" of any two values are 50% for velocity, 
and 45° for direction. In most cases, the angle values resulting 
for the correlated data are in agreement within 20 °, which is 
excellent. 

The 102 remaining events, unusable for the two indepen- 
dent velocity measurements, suffered from various difficulties. 
About 407o of the events were typified by very fast fluctuations 
in the insolation which resulted In difficulty identifying the 
event at corresponding stations . Another 40% of the events 
were generally well defined, but were missing at one of the sta- 
tions du'.e to shape, size, or position of the cloud over the field. 
Another 10% of the event losses were due to paper drive problems. 
Another 5% were lost due to problems with the electronics caused 
by electrical storms, and the rest were due to loss of event tim- 
ing and alignment clock marks due to power outages . 

The most common reason for disagreement between the two 
vector velocity calculations in the 93 ’ fully represented events 
was sign errors in the event timing. The events which result 
in these sign errors, do not appear to he correlated in cloud size 
(event duration) . They also do not correlate with days which 
presented difficult alignment or event Identification. The 
most likely cause of these errors is the random nonsymmetric 
shapes or peculiar angles of the leading and trailing edges with 
respect to the direction of motion. Of the 93 events, 62 events 
or 67 % show correlated vector velocities. It is our belief that 
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th-LS level of spatial correlation in the data is excellent con- 
sidering the complexity of the atmospheric effects involved in 
cloud motion and the simplicity of the mathematical model assumed 
for the calculations. In the following section, the results of 
these calculations with the correlated data are presented. 

3 . 3 Reduced Data - Velocity, Angle, and Size 

3.31 Correlated data 

An examination of the correlated data shows the findings 
regarding the local microclimate cloud patterns. These results 
are also compared with information available from the National 
Weather Service and the Institute of Atmospheric Physics at the 
University of Arizona. 

The correlated data shown in Figures 11-13 resulted in 
most probable cloud speeds between 12.8 and 16 Icm/hr (8 and 10 
mph) , a most probable direction of motion to the North-West or 
South-East, and a most probable size around 1.6 km (1 mile). 

Data published by the National Weather Service for summer Tucson 
weather, list average observed ground wind speeds during the 
months of June-August as approximately 8 mph. If we assume the 
upper level wind speeds are of the same order of magnitude, then 
we find excellent agreement for our data. Regarding the direction 
of motion of the clouds, it is a long recognized fact, confirmed 
by direct communication with the staff of the Institute of Atmos- 
pheric Physics, that summer storms in the Tucson Valley tend to 
move from the SE to the NW. This prevailing pattern is due to 
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Figure 11. Velocity distribution determined from 
spatially correlated data. 
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Figure 13, Distribution of calculated cloud sizes, 

D, (kilometers) from spatially correlated 
data . 
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the topological environment around Tucson, and the availability 
of water vapor from the Gulf of Mexico, The Helio detector field 
is located in the extreme Northeast corner of the Tucson Valley, 
within a few miles of the mountains on the East and the North. 
These mountains, with peaks higher than 2 700 m (9 000 ft), can 
alter the prevailing storm motions by significant amounts. The 
data shows a most probable direction of about 130° for the 
observed clouds in excellent agreement with the above pattern. 

The cloud size distribution, shown in Figure 14, shows 
a relatively strong probability for clouds of approximately 
1.6 km (1 mile) extent. This result agrees very well with 
qualitative visual observations of the cloud sizes. There 
appear to be many clouds of extent less than 0.S km (0.5 miles), 
but these do net appear in the data because of the time resolu- 
tion of the measurement system. Similarily, these observations 
indicate that individual clouds with sizes greater than 9_ lcm are 
rare. The specific numbers for velocity, direction and size 
quoted from the data are only approximate. The variance in the 
data due to the nature of the events being measured limits the 
accuracy of the determinations. 

3.32 Correlated plus uncorr e lated data 

To further amplify on the quality of the data as rep- 
resentative of the actual Tucson microclimate, we can examine the 
complete body of data, including the 134 uncorr datable events. 

If these data agree with the correlated data with respect to 
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cloud speed and direction, we can be confident of the quality 
of the data for use in the modeling. 

The speed and direction distributions for the complete 
set of data are shown in Figures 14 and 15. The most probable 
speed is still about 12.8 km/hr (8 mph) , but the most probable 
direction is closer to due West. Considering the latitude in 
the angle measurements and the possible direction alterations 
due to mountains, as discussed above, these results agr.eee with 
both the correlated data and the prevailing storm patterns in 
Tucson. 

3 . 4 Power Plant Site Analysis 

To analyze the microclimate for a particular power 
plant location there are several questions which must be answered. 
If we assume, for our model, that the insolation at the site is 
occulted by Individual clouds on days of partial cloudiness, 
then the questions that bear the most importance to the operation 
of a power plant are as listed in Table 2. The answers to 
these questions can be related to the possible design features 
of individual components of a collector system, or to the dynam- 
ics of the whole system. Let us use the central receiver photo- 
thermal process as an example. Questions 1-3 are relevant to 
the Initial placement of the system in an area with the least 
possible insolation Interference. Question 4 can give the 
designer an idea of the extent of dead time to be expected for 
each individual mirror in the concentrator field for partly 
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Figure 14. Velocity (speed) distribution for cloud motion 
over the field using all of the data. 
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Table 2. Microclimatology questions pertinent to power 
plant site analyses. 


ASSUMPTIONS: 1.) Each occultation is due to an individual 

cloud of approximately circular shape. 

2.) For power plant sizing, the system is 

assumed to have a 25% overall efficiency 
for a peak Input of 1000 w/m z , and the 
field dimension is adjusted to compensate 
for shading. 


QUESTIONS : 


1. ) How often do occultations occur. 

2. ) What is the size distribution for Indi- 

vidual occulting clouds, 

3. ) Is the temporal distribution for occult- 

ing events random or biased. 

4. ) How long will an elemental area in the 

field be effected. 

5 . ) What fraction of the field is affected 

by any single event. • 

6. ) How long will the whole field be effected 

by any single event. 
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cloudy days, and questions 5 and 6 relate to the magnitude and 

duration of differential inputs to the receiver caused by field 

’ ‘ ■ 

shading. The answers to these questions are presented below. 
The data and analysis results amply demonstrate the usefullness 
of these analyses. 

If we examine ’the data from our field of sensors in 
terms of the number of partly cloudy days vs the total number 
of days observed, we find that 29% of all Tucson days have 
cloud activity leading to possible short duration occultations. 
For this purpose, a partly cloudy day is assumed to be any day 
where cloud interference was observed for more than 15%, but 
less than 85% of the daylight hours. (See the discussion of 
the calendar day data as listed in Table 4) • If we examine 
the data for statistical bias, we find that there is a slightly 
greater probability that occultations will occur in the after- 
noon and evening hours . This is evident from looking at the 
calendar data at the end of this report. (On these records, 
the time of day reads from right to left prior to March IS, and 
left to right thereafter) . 

The size distribution for the events observed In Tucson 
was presented In the preceding section. The interpretation of 
this data with respect' to the effects on an elemental area are 
presented in Fig. 16. The data displayed here is the distri- 
bution of the ratio of cloud size to the velocity for each 
individual event. Therefore , this "dead time" represents the 
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Figure 16. Frequency vs the time any cloud is over 
an elemental area, AA, in the field. 


length of time that a particular cloud shades an elemental 
area of the field or an individual concentrator element. As 
the data show ; , the dead time is strongly peaked around 6 min- 
utes. The importance of this information for the central 
receiver concept, for example, is in the design of the individ- 
ual concentrator element tracking systems to minimize the loss 
of mirror alignment after the event clears . 

If we examine the fraction of the field area lost 
during each event, we find the results shown in Fig. 17. For 
large field systems where the size of the cloud is greater 
than the mean cloud sizes, the loss of receiver input from the 
field could cause differential heating problems. For the 
smaller systems, it is likely that the whole system would lose 
input due to the cloud, and that the system cooling would be 
a problem. 

The severity of both of these problems would depend 
on the length of the dead time associated with the event. These 
results — the answer to question 6 -- are presented in Fig. 18. 

As expected, the larger the field, the longer the occult at ion 

i * 

will cause problems. Data like displayed in Figs. 17 and 18, 
for any area intended for powerplant siting should he used in 
optimum size calculations and in the system design considerations. 
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Area of Cloud/ Area of Field 


Figure 17. Frequency vs the percentage of the field 
area that is occulted by an individual 
cloud. ■ . 




^down (hrs > 

Figure 18. Frequency vs down time for passage of 

cloud over whole field - largest dimension. 
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3.5 Fourier Representation 

We have found that the solar insolation days for Tucson 
fall into one of three categories. Approximately 30-35% of all 
days have cloudless skies or very little interference in the 
insolation due to clouds. Another 10-15% of all days have heavy 
cloudiness or more than 80% sky cover during the daylight hours. 
The remainder of the days have moderate cloud activity which may 
occur at different -times and alter the solar insolation for 15% 
to 75% of the day. The first type of day is relatively straight 
forward and simple to explain by the existing' insolation models'. 
These latter types of insolation days may be easily represented 
in terms of the frequency spectrum of the temporal changes in 
the insolation. A convenient mathematical model for this fre- 
quency representation is the Ttfell known Fourier representation 
of time dependent functions in the form 

ca 

<Kt) = %A o + l (A^cos (mrt/T) + B^sinCmrt/T) ) 


n 


where 


A = % j<$> (t)cos (mTt/T)<fit 


n 


-T 


B = Yp /<> (t)cos (mTt/T)<f>t. 
-T 


In our previous studies , we have found characteristic 
shapes for the data when expressed in this form. A plot of 
with respect to time" 1 , results in curves similar to those 
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shown in Figure 19 * These curves are representative of nearly 
all clear and partly cloudy days observed in Tucson at any 
particular site and any season of the year. Of course there are 
variances in the magnitude of the higher frequency components 
depending on -the rapidity of the fluctuations in the insolation. 

The Fourier expression may also be extended to show the 
temporal nature of -the cloud motions across a collector field. 

If we use the velocities calculated from the event analysis, 
in conjunction with a representative distance across the field, 
r.., one can arrive at a -model for the temporal-spatial insola- 

4 

tion measurements, or local microclimate of the form 

co X • • 

♦ <t.r) = %A q + £C^cos~(t + -f- + Y n ) ■ ■ 

In this form, the spatial characteristics of the insolation are 

represented by an additional phase term. 

To evaluate the applicability of this model to our data, 

we studied the ability to correlate the data between the five 

individual stations via the addition of various size phase terms 

to the Fourier representation of a particular day’s record. The 
* 

modeling was considered a success in those cases where the 
’’events” were reproduced from location to location with well 
preserved widths and depths, but with appropraite shifts in 
time. The time shifts were deemed excellent if velocities 
could be calculated and correlated as discussed above. 
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Figure 19. The Fourier Spectra for clear and partly 
cloudy days . 



The days which showed good Fourier spatial representa- 
tion ,were those days with relatively broad and well-defined 
oc.cultations ; by broad we mean events of duration greater than 
about 0.1 hr. It was noticed that data records which contained 
large amount s of high frequency structure, or many events with 
durations less than 0.1 hr length, were hard to correlate 
between station locations. These types of days usually were 
associated with a great deal of variance in both the breadth 
and depth of individual occultations between the five stations. 
This variance was not selective with respect to size of the 
event. The velocity correlations described above were also 
difficult for such days due to splitting of "events" or 
fading of an "event" between stations. In particular, occur- 
ances such as absence of an "event" on one or more of the 
records, or "splitting" of a broad event from one station to 
the next by a spike in the center of the dip, or the combining of 
two events into one broad dip, were common problems for days 
with appreciable high frequency activity. These types of 
phenomena are probably correlatable to particular types of 
cloud structure. It is most likely that the presence of many 
short duration events would correspond to light cirrus clouds, 
cumulus clouds with rough (wispy) eges, or other cloud struc- 
tures associated with extreme turbulence. The better events 
are most likely associated with well-defined and well-localized 
clouds . 
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The Fourier plus phase representation for each day was 
judged as good, good-fair, fair-poor, or poor depending on 
the quality of the fit to all stations. Of the 41 days 
analyzed, a good fit was found for only 15 days. If the fair 
to good days are included, the count Is Increased to 26 days. 

The fair to good days are those where the fit was applicable 
for more than half of the Identifiable events. The greatest 
fault was in the reproduction of the widths and depths of the 
events. The poor days were generally associated with large 
variences in both width and depth of the events or with missing 
events on one or more station record. The ability to relatively 
well account for the five stations, 26 out of 41 times using the 
phase shift model yields a 63% confidence in the model. We 
would, therefore, conclude that the technique is applicable to 
microclimate .predictions, but limited in its accuracy. 

A good estimate of the phase term can be caluculated 
by using the most probably magnitude of velocity as detemined 
from the velocity profile such as shoun In Figure 14. If one 
then references to a particular location in the field, the 
phase term can he calculated for any other point at a distance 


j r j away . 

By examining the Fourier representations of the days, 
one can type days according to the six classifications shown in 
Table 3. The applicability of the Fourier-Phase model for 
each type of day is shown by a "G” for good, a n G-F" for good 
to fair, or a "P" for poor. 
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Table 3. Day classifications and the results of Fourier plus 
phase modeling for each. 


DAY 

FOURIER ( 

TYPE 

DESCRIPTION 

+ 

NO. 


PHASE 

1 

Nearly clear with well defined occultations 
and very little high frequency activity. 

G 


Example: June 29, 1975 - Figure 27 


2 

Nearly clear with well defined occultations 
and some high frequency activity. 

Example: February 16, 1975 - Figure 23 

G-F 

3 

30-80% of day contains occultation "events*’ 
or activity. Much high frequency activity. 
No overcoat period. 

P 


Example: January 9, 1975 - Figure 22 


4 

30-807o of day contains occultation activity. 
Overcast for part of day. Little high fre- 
quency activity. 

Example: January 7, 1975 - Figure 22 

G 

5 

30-80 % of day contains occultation activity. 



Overcast for part of day. Much high frequency 
activity. 

P 


Example: January 27, 1975 - Figure 22 


6 

Overcast day with sharp peaks. 
Example: March 8, 1975 - Figure 24 

F-G 
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3 . 6 Calendar Representation 

The data obtained from the MBOR instruments can also 
be displayed in calendar format. This display shows, at a 
glance, the nature of the solar insolation during the months. 
The number of clear days, the trends through the weeks, and 
the types of cloudy days during the month, are all readily 
obtainable from even a cursory study of the calendar. 

Included in this report as Appendix 1, Figures 20-27 
are the calendar data for November 1974 through June 1975. 

After June, the Home Station recorder was shifted to the 2'7hr ' 
speed, making the calendar photography impractical. Inspection 
of these records yields some very interesting information. Let 
us assume the classification of the days as follows : 

•Clear - <15% of day's insolation 
interruped by clouds; 

Partially Cloudy - >157,, but < 8 5% 
of day interrupted; 

Overcast - >857. of day interrupted. 

Application of this scheme to the five months data results in 
the breakdowns shown in Table 4- The data shows a rather 
remarkable constancy in the fractions of the type of insolation 
days observed during the three seasons. 
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Table 4. Diyision of Calendar-Day data into three classifica- 
tions showing similarity of the eight Winter /Spring/ 
Summer months . 


NO. DAYS 
OF DATA 

. MONTH 

‘ NO. OF DAYS 

o.c. 

NO. OF DAYS 
PARTIAL 

NO. OF DAYS 
CLEAR 

30 

Nov. 

1974 

2+7% 

11+37% 

17+57% 

30 

Dec. 

1974 

3+10 % 

5+17% 

22+73% 

31 

Jan. 

1975 

3+10% 

8+25% 

20+65% 

26 

■Feb. 

1975 

1+4% 

8+31% 

17+65% 

30 

Mar. 

1975 

2+7% 

12+40% 

16+53% 

26 

Apr. 

1975 

5+17% 

9+30% 

12+40% 

31 

May 

1975 

1+3% 

10+32% 

20+65% 

26 

Jun. 

1975 

0+0% 

4+13% 

22+73% 
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3. 61 Comparisons with. MTS data 

We may also present the data in Table 4 in graphical 
form as shown by the solid lines in Figs. 28 and 29. The data 
can then be compared to the pertinent National Weather Service 
information for the same months. The two sets of NSW data 
which should have close correlation to our "days type" analysis 
are the percent of cloud cover, and the percent possible sun- 
shine. The observations of percent cloud cover are of a quali- 
tative nature because they depend on observational judgements of 
individuals regarding the fraction of the slcy. that contains 
clouds. It is our experience that the proximity of mountains 
and the "bowl" shape of the Tucson Valley yield many days where 
there is significant amount of cloud cover which only minimally 
obscures the valley floor. Much of the cloudiness is well 
localized over the mountains in these cases. This may be the 
reason for the relatively poor agreement in the percentages 
displayed in Fig. 28. In this figure,, the lower lines are the 
percentage of days fo.r each month that correspond to "clear" 
skys. Clear skys are defined as days with ’less than 15% activity 
or less than 1 tenth cloud cover. The area between the two lines 
represents the partly cloudy percentage given by 15-85% activity 
and 1-9 tenths cloud cover. 

Of perhaps more significance is the correlation between 
the "percent of possible sunshine 1 ’' and our analog data. The % 
sunshine readings are calculated from the ratio of the total 
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HELIC' Data (Solid Line ) 
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Partly Cloudy => 15%< Activitij <85% 
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NWS % Cloud Cover (Dashed Line) 

Clear => 0-10% 

’Partly Cloudy => 20<% Cloud Cover <90% 

Overcast => 100 % 


Figure 28. Percentages of days elassifed as clear, partly 
cloudy, or overcast for the months of November 
1974 through June 1975. 
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HFLIO Data (Solid Line) 

Clear . => <15% Activity 

Partly Cloudy —> 15%< Activity <85% 

Overcast => >85% Activity 

NWS % Possible Sunshine (Dashed Line ) 

Clear => >85% 

Partly Cloudy => 15<% Sunshine <85% 

Overcast => <15% 


Figure 29. Percentages of days classified as clear, partly 
cloudy, or overcast for the months of November 
1974 through June 1975. 




number of hours when the direct insolation is greater than zero 
to the number of hours between sunrise and sunset. There are 
many times when occulting clouds to not totally obscure the sun, 
i.e. the direct flux does not go to zero during the occultation 
Examples of this are apparent in the calendar records for Novem 
ber In Fig. 20. November 13 and 14 are excellent examples. 
These difficulties could easily account for the differences 
shown in Fig, 29.. The interesting result of our measurements 
is that the partial occultations are evident in these records. 
The 7c possible sunshine- data can be significantly misleading 
as to the total direct insolation available at a site. It 
should be emphasized that the device used to record the "possi- 
ble" sunshine record is merely a switch which turns on or off 
as the insolation exceeds or drops below a threshhold value. 

The device is usually adjusted to "turn-on" when the direct 
flux is just enough to yield a definite shadow for an object. 
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SECTION 4,0 CONCLUSIONS 

The analysis of the data as presented above yields much 
valuable information for the designer of large scale solar 
absorber systems. The data obtained from this study can be 
used to develop the most efficient system configuration and/or 
controls for any system for any particular site. The micro- 
climate statistics available from these measurements are a good 
representation of the insolation dynamics expected for a large 
field array. This data can be statisically analysed to allow 
design input of the probability of insolation Interference and 
the dead times expected for individual elements of the field 
or the field as a whole. These data, in conjunction with the 
area losses and imbala:ices due to occultations, are essential 
to the achievement of a best cost effective design for large 
solar collector systems. 

Comparisons of the Hello data records to the NSW data, 
results in qualitative correlation. The differences in the 
two data sources seem to be attributable to the loss of infor- 
mation content inherent to the methods of measurement employed 
by the NWS or to the separations of several miles between the 
NWS station and the Helio network. 
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SECTION 5,0 RECOMMENDATIONS FOR CONTINUED RESEARCH 


The statistical analyses presented here have demonstrated 
the usefulness of the microclimatological modeling of power plant 
sites. The data base presented is, however, only minimal, and 
should be expanded and refined. It is our belief from visual 
observations that a great deal of the cloud occultations had 
periods less than the temporal resolution of the system. It Is, 
therefore, recommended that this resolution be improved. This 
is easily accomplished by increasing the recorder drive speeds 
to "spread out" the data, and the improvement- of the clock 
timing system used to correlate the five individual records. 

In addition, there are other data correlations which 
may prove very useful. In particular, the use of correlated 
photography of the occulting clouds and their type may be 
useful in relating the field temporal statistics to cloud type 
$nd cloud cover statistics available from the Weather Service. 

The data analysis should also be expanded to yield comparative 
models for both focusing and non-focusing systems in which we 
may be able to describe the magnitude of the insolation changes 
at any point in the field, A^CR^), in terms of a probability 
distribution. These distributions could be calculated for 
both the direct and the total insolation for applications to 
concentrating or planar systems. This modeling would present 
data on the fractional loss of input energy and would be 
useful in specifying the type of system best suited to a 
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particular site. The continuation of this worh is the subject 
of a proposal presently being drafted. 
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APPENDIX 1 


Figures 20 through 27 : 

Calendar display of the analog data for the months 
of November 1974 (Fig. 20), December 1974 (Fig. 21), January 
1973 (Fig. 22), February 1975 (Fig. 23), March 1975 (Fig. 24), 
April 1975 (Fig. 25), May 1975 (Fig. 26), and June 1975 (Fig. 27). 
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